TWr METHOD FOR AN ENTROPY FILTER ECHO PROCESS.NG IN 
' TIME-OF-FLIGHT OR LEVEL MEASUREMENT SYSTEMS 



EIF1 n OF TH P INVENTION 

10001] The present invention relates to signal processing, and more 
particularly to a method for an entropy filter for echo processing in level 
measurement or time of flight ranging systems. 



BACKGROUND nF ruF - INVENTION, 

[0002] Pulse-echo acoustic ranging systems, also known as time-of-fl.ght 
ranging systems, are commonly used in level measurement applicat.ons. 
Pulse-echo acoustic ranging systems determine the distance to a reflector 
(i e reflective surface) by measuring how long after transmission of a burst 
of energy pulses the echoes or reflected pulses are received. Such systems 
typically use ultrasonic pulses, pulsed radar or microwave s,gnals. 



[0003] Pulse-echo acoustic ranging systems generally include a 
transducer and a signal processor. The transducer serves the dual role of 
transmitting the energy pulses and receiving the reflected energy pulses or 
echoes. An echo profile is generated from the received energy pulses. Echo 
pulses are identified in the echo profile by the signal processor, and the 
distance or range of the object is calculated based on the transmit times of 
the transmitted energy pulses and the receive echo pulses. 



[0004] To provide accurate level measurements, the echo pulses must be 
precisely detected and processed. One approach involves transm.tt.ng 
multiple transmit pulses, i.e. "shots", and averaging the echo pos,t,ons 
associated with the multiple transmit pulses. For noisy environments, for 
example a waste water tank, the number of shots is increased for the 
average in an effort to reduce the effects of noise. The aim of this approach 
is to dampen or smooth random deviations in the echo pulse position in an 
echo profile. However, the technique does not prevent noisy profiles from 
becoming part of the average. In a noisy environment, such as filling a waste 
water tank, the level measurements can be very imprecise because the 
average echo-position is based on a large number noisy echoes. 



[00051 The commonly used technique for finding echoes in an echo profile 
involves generating a time varying threshold or TVT curve. The TVT curve 
provides a baseline or line on the echo profile which is above the noise level 
in the echo profile. Valid echoes appear above the TVT curve. Various 
algorithms and techniques are known in the art for determining the no.se 
floor and generating the TVT curve. 



[0006] Two solutions for improving level measurements for echo-based 
systems in the presence of noise comprise: (1) increasing the number of 
shots taken and averaging the results; and (2) measuring the noise floor and 
using the average noise level to generate a TVT curve to assist in the 
interpretation of the echo profile curve. 



[00071 A typical echo profile indicated by reference 1 00 is shown in Fig. 1 
together with a TVT curve indicated by reference 120. The first portion of the 
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filter echo profiles that do not have a minimum amount of bit-i 
information, i. . a low degree of randomness. 



[0011] In a first aspect, the present invention provides a method for 
generating an echo signal in a time-of-flight ranging system, the method 
comprises the steps of: transmitting one or more bursts of energy towards a 
reflective surface; receiving reflected pulses from the reflective surface, and 
converting the reflected pulses into an echo signal, the echo signal includes 
one or more potential echo pulses; applying an entropy filter to the echo 
signal, the entropy filter includes determining whether information contained 
in the echo signal corresponds to a valid echo pulse or is substantially noise; 
distinguishing those of the echo pulses determined as comprising noise in 
the echo signal. 



[0012) In another aspect, the present invention , provides a level 
measurement device for measuring a distance to a material having a 
surface, the level measurement device comprises: a transducer for emitting 
energy pulses and detecting energy pulses reflected by the surface of the 
material; a controller having a receiver and a transmitter; the transducer 
includes an input port operatively coupled to the transmitter and is 
responsive to the transmitter for emitting the energy pulses, and the 
transducer Includes an output port operatively coupled to the receiver for 
outputting reflected energy pulses coupled by the transducer; the receiver 
includes a converter for converting the reflected energy pulses into signals; 
the controller includes a program component for generating an echo profile 
based on the signals, the controller includes another program component for 
applying an entropy based filter to the echo profile. 



W0131 In a yet further aspect, the present invention provides a time of 
flight ranging system comprising: a transducer for emitting energy puis s 
and detecting reflected energy pulses; a controller having a receiver and a 
transmitter; the transducer includes an input port operatively coupled to the 
transmitter and is responsive to the transmitter for emitting the energy 
pulses, and the transducer includes an output port operatively coupled to the 
receiver for outputting reflected energy pulses coupled by the transducer; 
the receiver includes a converter for converting the reflected energy pulses 
into signals; the controller includes a program component for generatmg an 
echo profile based on the signals, the echo profile comprises potenual 
echoes and potential noise; the controller includes another program 
component for applying an entropy Alter to the potential echoes to 
differentiate the potential noise from said potential echoes. 



,0014] Other aspects and features of the present invention will become 
apparent to those ordinarily skiiled in the art upon review of the following 
description of specific embodiments of the invention in conjunction w,th the 
accompanying drawings. 



BRIEF DESCRIPTION OF TH E DRAWINGS 

[00151 Reference is next made to the accompanying drawings which 
show, by way of example, embodiments of the present invention and m 



which: 
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[0016] Fig 1 is a graphic representation of an echo signal waveform, a 
noise floor level waveform and a noise signal waveform for a time-of-fl.ght 
ranging system; 

10017] Fig. 2(a) is an echo profile waveform with echo pulses identified in 
accordance with the present invention; 



[0018] Fig. 2(b) is an exemplary echo profile with a defined or valid echo 
pulse; 

[0019] Fig. 2(c) is an exemplary echo profile comprising mostly noise with 
no clearly defined echo pulse; and 



[0020] Fig. 3 is a flow chart showing the processing steps for an entropy 
filter for processing an echo profile in accordance with the present invention. 

n ft All FD DESCRIPTION OF EM BODIMENTS OF THE INVENTION, 

[0021] Reference is first made to Fig. 2(a) which shows an echo profile 
200 generated in a time-of-flight ranging or level measurement system and 
suitable for processing in accordance with the present invention. 

[0022] The echo profile 200 is generated in the operation of the level 
measurement device (or the time-of-flight ranging system). In known manner, 



the level measurement device includes a transducer (e.g. ultrasonic, 
microwave or radar), a controller or signal processor unit, an analog-to- 
digital (A/D) converter, a transmitter, a receiver, and a power supply unit. 
The transducer emits a transmit pulse or energy burst directed at a surface 
to be measured. The surface reflects the transmit energy burst and the 
reflected energy pulses are coupled by the transducer and converted into 
electrical signals. The electrical signals are applied to the receiver and 
sampled and digitized by the A/D converter. The signal processor, for 
example a microprocessor operating under firmware control, takes the 
digitized output and generates the echo profile 200 having a form as shown 
in Fig. 2 The echo profile 200 is characterized by one or more valid echoes 

202, Indicated individually by references 202a, 202b, 202c The echoes 

202 correspond to reflected energy pulses. The echo profile 200 also 
includes a ring-down portion 204 which corresponds to the 'ringing down' of 
the transducer after transmission of the energy burst or pulses. It is normally 
not possible to detect receive echoes during the ring-down interval. As also 
shown in Fig. 2(a), the echo profile 200 includes a multitude of pulses or 
bumps 206. indicated individually by references 206a. 206b, 206c, 206d, 
206e, 206f,... which are identified as noise. As will be described in more 
detail, the noise pulses or elements 206 in the echo profile 200 may be 
distinguished by utilizing the entropy filter in accordance with the invention, 
or by applying a TVT analysis in conjunction with the entropy filter according 
to the present invention. 



[0023] In known manner, the controller unit executes an algorithm (e.g. a 
process or function implemented in firmware) which uses the echo profile 
200 to calculate the range, i.e. the distance to the reflective surface, from 
the time it takes for the reflected energy pulse to travel from the reflective 
surface to the transducer. From this calculation, the distance to the surface 
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of the liquid (i.e. the materia! contained in the vessel) and thereby the l< 
of the liquid is determined. The controller, e.g. microprocessor 
microcontroller, is suitably programmed to perform these operations as 
be within the understanding of those skilled in the art. 



[0024] As will now be described in more detail, the subject invention is 
directed to a method or process for an entropy filter for analyzing an echo 
profile for echoes or pulses. The method and processing steps as described 
below may be embodied in the controller as a program component or 
firmware. 



[0025] In operation the entropy filter distinguishes or rejects echo profiles, 
or "shots', that do not contain a minimum threshold amount of information, to 
define a valid echo pulse, on the basis that the more disordered the bit 
stream in the echo profile, the less reliable the echoes (i.e. shots) for 
measurement, regardless of the magnitude or threshold level of the echoes. 
An exemplary echo profile 220 having a defined echo pulse 222 is depicted 
in the Fig. 2(b). In Fig. 2(b), a TVT curve 224 is also included on the plot for 
the echo profile 220. Another exemplary echo profile 230 is depicted in Fig. 
2(c). The plot in Fig. 2(c) also includes a TVT curve 234. As compared to the 
echo profile 220 of Fig. 2(b), the echo profile 230 of Fig. 2(c) is a much more 
noisier profile with a majority of pulses 232, indicated Individually as 232a, 
232b, 232c, 232d, 232e and 232f, falling below or around the TVT curve 234. 
As a result it Is very difficult to ascertain a clearly defined echo pulse for the 
echo profile 230 of Fig. 2(c). According to this aspect of the Invention, 
application of the entropy filter to the echo profile 230 of Fig. 2(c) determines 
more random information which manifests itself in irregular contours, e.g. the 
noise pulses 232. The high degree of randomness is indicative of noise as 
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compared to valid echo pulses. For the echo profile 220 of Fig. 2(b). 
application of an entropy filter in accordance with the present invention 
determines higher information content (i.e. a low degree of randomness) and 
the echo pulse 222 is identifiable. 



[0026] in accordance with this aspect, the entropy of an echo profile (e.g. 
the echo profile 220 shown in Fig. 2(b)) is determined by match.ng a 
reference set of short bit-patterns with selected portions of the bit-stream 
comprising the echo profile to provide an information content estimate for 
the echo profile. The portions of the bit-stream for the echo profile may be 
randomly selected. 



[0027] One technique for analyzing a random portion of the echo profile 
involves applying a Monte Carlo method. Using the Monte Carlo method, if 
each reference pattern occurs the same number of times in the selected 
portion then the selected portion, e.g. echo profile, has the minimum amount 
of information possible. In addition, the more that some reference patterns 
appear in the bit-stream, the more information content and therefore less 
entropy i.e. randomness, the echo profile is estimated to have. Th.s 
increases the likelihood or probability that the echo profile comprises an 
accurate or valid echo position. 



[0028] According to this aspect, the echo profile is provided as a b.t- 
stream comprising a binary sequence represented by X and containing K 
possible patterns or symbols, where each pattern x, occurs with probably 
p.. The entropy of the bit-stream, i.e. echo profile, is defined by a function 
H(x) as follows in equation (1): 
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K 

HM = - 5LP*log 2 p* (D 
k-1 



[0029] If a symbol or pattern occurs with probability 1, for example, the 
bit-stream for the echo profile comprises I's, then the entropy is at a 
minimum, i.e. H(x) - 0. If there are K possible symbols or patterns, and each 
of the symbols occurs with equal probability, then the entropy is at a 
maximum, i.e. H(x) = log^K and the echo profile comprises noise. 

[0030] For example, for an echo profile X, the entropy is determined by 
matching 14 symbols, i.e. K=14 comprising short binary patterns, x k , with 
(e.g. randomly selected) portions of the echo profile. In this example, the 
short binary patterns, x k , comprise 0, 1, 00, 01, 10, 11, 001, 010, 100, 110, 
011, 101, 000, 111. The echo profile X comprises a bit-stream representing 
echo profile magnitudes which are binned with respect to a radix and 
precision. 



[0031] According to another aspect, each echo profile, for example, a 
plurality of echo profiles generated from a batch of 'shots" for a 
measurement, may be ranked according to its entropy. The rankings may 
then be used by processes in the application software for the level 
measurement system to select the most accurate echo positions for 
determining an average echo position and corresponding level 
measurement. 
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[0032] Reference is next made to Fig. 3 which shows in flow-chart form 
an entropy filter process 300 in accordance with the present invention. As 
indicated in block 302, the first step in the entropy filter process 300 
comprises selecting a portion of the echo profile (in the alternative the 
complete echo profile may be considered). The selected portion of the echo 
profile is represented digitally, i.e. as a bit wise data stream, with digital 
values representative of magnitudes in the echo profile. Next in block 304, 
one or more bit patterns are selected. The bit patterns may be generated or 
read from a look-up table stored in memory. Next in block 306, the bit 
pattern(s) are applied to the bit wise representation of the echo profile, and 
the number of occurrences of the bit pattern in the echo profile is 
determined- As described above, the Monte Carlo method or technique may 
be utilized as described above. The number (or probability) of occurrences of 
the bit pattern(s) indicates the entropy of the echo profile. A high number of 
occurrences corresponds to a low entropy and indicates an echo profile with 
low noise content and therefore a valid echo pulse. A low number of 
occurrences corresponds to a high entropy which is indicative of noise. 
Based on the number of occurrences (decision block 308), the echo pulse 
profile is characterized as an echo profile having a defined echo pulse (block 
310) or as noise (block 312). If more than one echo profile is being 
considered, for example, to generate an average echo profile using 
application of the entropy filter to a number of echo profiles (i.e. shots), then 
the process is repeated (decision block 314), otherwise the process returns 
(step 316) to the calling function or process being executed by the controller. 
It will be appreciated that the processing steps illustrated in Fig. 3 are 
exemplary, and as such may be implemented in a different fashion or order. 
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[00331 An exemplary implementation in the c programming language is 
provided as follows. The specific implementation details of the code will be 
within the understanding of those skilled in the art. 



Pseudo code 

© Siemens Milltronies Process Instruments Inc. 2002-2003 

float calculateShotEntropy(void) 

{ 

r randomly select a section of profile data */ 
profileOffset = (size j)(profileSize * rand()); 
sectionStart = profileOffset; 
secttonEnd = profileOffset + sectionSize; 

/- convert profile section to binary */ 

X = convertToBinary(profileData, sectionStart, sectionEnd); 

I' calculate entropy for the section V 
entropy = 0.0; 

for (k = 0; k <; numberOfSymbols; k+*) 
{ 

/« sum number of occurrences of this symbol in section */ 
numOccurences = 0; 

for (bitOffset = 0; bitOffset < sectionSizelnBits; bitOffset++) 
{ 

if (matches(symbol[k], X[bitOffset]) 
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{ 

numOccurrences++; 
} 

} 



, probab l.lty for «. sy.bo, occurring is number of occurrences divided b y number 
of possible occurrences */ 
ptkl = numOccurrences / (sectionSizelnBits - 1); 
r Add entropy for this symbol to running total for entropy V 
entropy +* -(p[k]*«og(p[kl,2)); 
} 

return(entropy); 

} 

,0034, ,« will be appreciated «ha. .he entropy shot filter processing in 
accidence with the present invention complements and/or improves the 
S trom avera g ,n 9 batches o, shot data. Averaging is an a«emp. 
sm „oth out random perturbations in the echo profiie, however^ , p,c t ee 

averaging tends to only remove noise that has a 

distribution. The entropy filtering technique according to the present 
vlon detects noise with arbitrary probability distributions acoordm to 

he symmetries of the selected set of reference bit patterns. Furthermore 

Z enlpy fitter processing may be utilized to complement existing no.se 
n Lnioues. for example, caicuiation o, a noise hoor for generating 



a TVT. 



[00351 The present invention may be embodied in othe spec hc _hm» 
lo « departing from the spirit or essentia, characteristics thereof Cert 
adaptations and modifications of the invenhon will be obvious to those skilled 
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in the art. Therefore, the above discussed embodiments are considered to 
be illustrative and not restrictive, the scope of the invention being indicated 
by the appended claims rather than the foregoing description, and all 
changes which come within the meaning and range of equivalency of the 
claims are therefore intended to be embraced therein. 



